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Introduction 

As the era of a permanently manned Space Sta- 
tion begins to unfold, the necessity for improved ra- 
diation protection methods for the astronauts who 
will man it increases in importance. For the first 
time in the history of the manned space program, 
the combination of long-duration missions and career 
astronauts will cause the high-energy heavy-ion com- 
ponent of galactic cosmic rays to assume major radio- 
biological significance. To properly assess spacecraft 
and personal shielding requirements and to prop- 
erly determine radiation exposures of critical organs, 
methods for estimating the radiation fields produced 
by relativistic heavy nuclei incident upon thick layers 
of absorbers of different composition are required. To 
meet these requirements, an extensive program for 
the development of theoretical methods for predict- 
ing relativistic heavy-ion interaction and transport is 
in progress at the Langley Research Center. Concur- 
rently, experimental verifications of these models are 
obtained by researchers at Lawrence Berkeley Labo- 
ratory using state-of-the-art equipment (ref. 1) and 
techniques (refs. 2 and 3). 

In previous work (refs. 4 to 8), a simple theory 
of heavy-ion fragmentation based upon a two-stage 
abrasion-ablation model has been developed for use 
in estimating the fragmentation parameters required 
for input into an accurate heavy-ion transport theory 
under concurrent development (refs. 3 and 9 to 11). 
Since the abrasion formalism calculates only the cross 
section for removal of a given number of abraded 
nucleons without specifying their identities as either 
neutrons or protons, a prescription for calculating the 
charge dispersions of the prefragments (the projectile 
pieces remaining after abrasion) is needed in order to 
calculate final isotopic and/or elemental production 
cross sections from the ablation process. Previously 
(refs. 6 to 8), two methods of determining prefrag- 
ment charge dispersions have been utilized in the 
Langley fragmentation model: (1) the hypergeomet- 
ric distribution, which treats the neutron and pro- 
ton distributions as completely uncorrelated (ref. 12), 
and (2) the giant dipole resonance (GDR) distribu- 
tion based upon the zero-point vibrations of the giant 
dipole resonance of the projectile nucleus (ref. 13). 
Both methods were used to calculate element produc- 
tion cross sections for 1.88 GeV/nucleon iron nuclei 
fragmenting in carbon, silver, and lead targets, and 
were also used to estimate isotope production cross 
sections for 213 MeV/nucleon argon nuclei fragment- 
ing in a carbon target. When estimating element 
production cross sections (ref. 6), better agreement 
between theory and experiment was obtained for the 
GDR method, whereas the hypergeometric distribu- 


tion yielded the best predictions for the isotope pro- 
duction cross sections (ref. 7). These later results, 
from reference 7, are displayed in figure 1, where it is 
apparent that the GDR calculations are too narrow 
in width. In this work, attempts to broaden the dis- 
tribution widths of these cross sections are made by 
incorporating a physically correct charge dispersion 
model based upon the Langley quantum mechani- 
cal abrasion theory (ref. 5) rather than the geomet- 
ric one used by Morrissey and coworkers (ref. 13). 
As will become apparent, the replacement of the 
geometric abrasion charge dispersion model by a 
quantum mechanical one yields only slight improve- 
ment in the agreement between theory and experi- 
ment for production cross sections of some isotopes 
produced in collisions of 213 MeV/nucleon 40 Ar with 
carbon targets. 

Symbols 


Af 

prefragment mass number 

Ap 

projectile nuclear mass number 

At 

target nuclear mass number 

B(e) 

average slope parameter of nucleon- 
nucleon scattering amplitude, fm 2 

b 

projectile impact parameter, fm 

C{ y) 

Pauli correlation function 

e 

two-nucleon kinetic energy in their 
center of mass frame, GeV 

GDR 

giant dipole resonance 

J 

droplet model coefficient (25.76 MeV) 

m 

number of abraded nucleons 

rh 

mean number of abraded nucleons 

N 

total number of nuclear neutrons 


renormalization coefficient for jth 
prefragment 

n 

number of abraded neutrons 

P(b) 

probability of not removing a single 
nucleon by abrasion 

Pr{Zi,Aj) 

prefragment charge-to-mass dispersion 
function 

Q 

droplet model coefficient (11.9 MeV) 

u 

defined in equation (12) 

Ax 

defined in equation (14) 

y 

two- nucleon relative position, fm 

z 

total number of nuclear protons 


Zp number of fragment protons 

Zp total number of projectile protons 

z number of abraded protons 

z 0 longitudinal position of projectile 

center of mass, fm 

ol z defined in equation (11) 

£r collection of constituent relative 

coordinates for target, fm 

pp projectile nuclear density, fm“ 3 

Pt target nuclear density, fm -3 

o-(e) average nucleon-nucleon total cross 

section, mb 

v*hT{Z)A) cross section for production of nucleus 
of type (Z, A ) by abrasion, mb 

&F fragmentation cross section, mb 

&m cross section for abrading m nucleons, 

mb 

Theoretical Development 


P(b) =exp j- A t cr(e) [2 tt B{e)\ 3 / 2 j dz Q 

x J d3 C T Pt($ t ) J d 3 y pp(b + z 0 +y + £ t ) 

(4) 

Methods for determining the appropriate nuclear 
distributions pi(i = P, T), constituent-averaged 
nucleon-nucleon cross sections cr(e), nucleon-nucleon 
scattering slope parameter B(e), and Pauli correla- 
tion function C{ y) are described in reference 6 and 
references eited therein. 

Hypergeometric Charge-to-Mass Dispersion 

The hypergeometric charge-to-mass dispersion 
function is based on the assumption that there is no 
correlation at all between neutron and proton distri- 
butions. Therefore, unphysical results such as abrad- 
ing all neutrons or protons from a nucleus while leav- 
ing the remaining fragment intact are possible. If z of 
the original Z projectile nucleus protons are abraded 
along with n of the original N projectile neutrons, 
then the hypergeometric distribution yields 


In order to ultimately calculate nuclear fragmen- 
tation cross sections, the cross section for forming 
a particular prefragment of mass Aj and charge Zj 
must be specified (refs. 6 and 7). In terms of the cross 
section for abrading m nucleons, < 7 m , the prefragment 
formation cross section is written as 

where 

^abr {Zj )Aj) = Pr {Zj , Aj ) a m (1) 

where Pr(Zj,Aj) is the prefragment charge-to-mass , 

dispersion function (to be specified) and is given and 
by (ref. 6) 


CTm = (i P ) / d2h t 1 " P W AF ( 2 ) 

where P(b) is the probability as a function of impact 
parameter of not removing a nucleon in the collision, 
and Afi the residual mass (prefragment) number, is 


with 


and 



( 5 ) 


A P = N + Z 

(6) 

m = n + z 

( 7 ) 

Zj — Z — z 

(8) 

Aj = Ap — m 

( 9 ) 


Ap = Ap — m ( 3 ) 

Within the context of eikonal scattering theory, an 
optical model potential approximation to the exact 
nucleus-nucleus multiple-scattering series, which in- 
cludes Pauli correlation effects, yields (ref. 6) 


As an alternative to the hypergeometric distribu- 
tion, Morrissey et al. (ref. 13) proposed a charge- 
to-mass dispersion distribution based upon the zero- 
point vibrations of the giant dipole resonance of 
the projectile nucleus. In this model, equation (5) 
becomes 


2 


Pr ( Zj,Aj ) = Nj (2ttc«?) 1 ^ 2 exp | 


~ [Zj-AjjZ/Ap)]' 
2a? 


( 10 ) 


where the variance (dispersion) is 




with 




(1 + u)- 3 / 4 

( 11 ) 


3J 

u = m 

Q(A P ) 1/3 


( 12 ) 


In the droplet model of the nucleus, the coefficients 
J and Q have the nominal values of 25.76 MeV and 
11.9 MeV, respectively (ref. 13). The rate of change 
of the number of nucleons removed as a function of 
impact parameter dm/db is calculated numerically by 
using the geometric abrasion model of reference 13. 
The normalization factor Nj insures that for a given 
value of Aj : the discrete sum over all allowed values of 
Zj yields unity for the dispersion probabilities. This 
overall normalization is not included in the original 
model of reference 3 but was added by us to insure 
probability conservation. 

Quantum Mechanical GDR Dispersion 

In an attempt to improve the GDR dispersion 
model, the method was generalized to a fully quan- 
tum mechanical treatment by replacing the geomet- 
ric abrasion model of reference 13 with the Langley 
abrasion model of references 5 and 6. Recall that in 
the geometric formalism, the variance was given by 
equation (11), which can be rewritten as 

— (£)(£)** < i3 > 

where from equation (11) we find 

/ n \ 1/2 

Az = 2.619 (— J (l + u)“ 3/4 (14) 

In references 6, 7, and 13, dm/db was calculated nu- 
merically by using a geometric abrasion model. In 
order to calculate this quantity quantum mechani- 
cally, we note that P( b) is the probability of not 
removing a nucleon by abrasion. Hence 1 — P( b) 
is the probability that a nucleon is removed. Clearly 
the mean number of nucleons removed as a function 
of impact parameter is then 

m = A P [l-P{b)\ (15) 


Therefore, the rate of change of the number of nu- 
cleons removed as a function of impact parameter 
becomes 


dm _ A dP( b) 
~db~~ P ~db~ 


(16) 


where P( b) is obtained from equation (4). In the 
actual calculation, values from equation (16) are ob- 
tained numerically. These values are then substi- 
tuted into equation (13) to determine the dispersion. 
The latter is then utilized in equation (10) to obtain 
values for the quantum mechanical GDR charge-to- 
mass dispersion model. 


Results 

To illustrate the results of these calculation meth- 
ods, cross sections for the production of sulfur, phos- 
phorus, silicon, and aluminum isotopes from the frag- 
mentation of 213 MeV/nucleon 40 Ar projectile nuclei 
by carbon targets are listed in table I. Displayed are 
theoretical predictions obtained from the fragmenta- 
tion model of references 6 and 7 for a hypergeometric 
distribution (eq. (5)) and for the GDR distribution 
(eq. (10)) evaluated for both' the geometric and quan- 
tum mechanical models. Also listed are experimental 
values estimated from the results of reference 14. As 
is apparent, only minor improvement in the agree- 
ment between theory and experiment is obtained for 
the quantum mechanical GDR model when compared 
with the geometric GDR model. This is clearly indi- 
cated in table II, where the x 2 values (ref. 15) for indi- 
vidual element production are listed for each charge- 
to-mass dispersion method. From tables I and II, it 
is clear that the hypergeometric distribution yields 
much better agreement with experiment than either 
of the GDR methods. The rather large x 2 value for 
aluminum production from the hypergeometric dis- 
tribution comes mainly from the relatively insignifi- 
cant 32 A1 isotope. If this isotope is excluded, the x 2 
value reduces to 8.4 for the hypergeometric distribu- 
tion, to 17.3 for the geometric GDR distribution, and 
to 16.7 for the quantum mechanical GDR distribu- 
tion. 

The results of this work confirm the findings of 
reference 12, in which it was noted that the hypergeo- 
metric dispersion relation yielded better agreement 
with experiment for predicted isotopic production 
cross sections than those obtained with a “complete 
correlation” function given by (ref. 12) 

Pr(Z j ,A j ) = ^m (17) 

This distribution, which would be correct if the nu- 
cleus were treated as a two-component crystal, is 
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not used in this work. Although not presented in 
this work, similar conclusions are obtained for the 
fragmentation results of 16 0 by Be and Pb targets 
with the charge-to-mass dispersion models described 
herein. Hence, future improvements to the Langley 
abrasion-ablation model will require improved meth- 
ods for estimating the prefragment excitation energy 
as well as the incorporation of final-state interactions 
into the theory (ref. 16). 

Concluding Remarks 

Methods for describing the charge-to-mass dis- 
persion distributions of projectile prefragments are 
described and used to determine individual isotope 
production cross sections for various elements pro- 
duced in the fragmentation of relativistic argon nu- 
clei by carbon targets. Although slight improvements 
in predicted cross sections are obtained for the quan- 
tum mechanical giant dipole resonance (GDR) dis- 
tribution when compared with the predictions of the 
geometric GDR model, the closest agreement be- 
tween theory and experiment continues to be ob- 
tained with the simple hypergeometric distribution, 
which treats the nucleons in the nucleus as com- 
pletely uncorrelated. 

NASA Langley Research Center 
Hampton, VA 23665 
December 13, 1984 
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TABLE I. ISOTOPE PRODUCTION CROSS SECTIONS FOR THE REACTION 40 Ar + 12 C — A Z + X 

[Incident kinetic energy is 213 MeV/nucleon] 



Isotope production cross sections, m 

b 

Isotope 

Hypergeometric 

Geometric GDR 

Quantum mechanical 

Experiment 

produced 

model (ref. 7) 

model (ref. 7) 

GDR model 

(ref. 14) 

37g 

7 

0 

0 

10 ± 5 

36 S 

31 

6 

6 

19 ± 8 

3S S 

34 

27 

25 

27 ± 10 

34 g 

60 

81 

81 

61 ± 20 

33 S 

35 

44 

43 

29 ± 10 

32 g 

10 

6 

6 

9 ± 4 

31 S 

.2 

.03 

.04 

1 ± 0.5 

36p 

0.3 

0 

0 


35p 

9 

1 

1 


34p 

15 

2 

2 

11 ± 4 

33p 

27 

25 

25 

25 ± 10 

32p 

30 

40 

40 

41 ± 10 

31 p 

16 

22 

22 

31 ± 9 

30p 

12 

3 

12 

9 ± 3 

29p 

0 

0 

0 


34 Si 

0.3 

0 ! 

0 

0.3 ± 0.15 

33 Si 

2.6 

0 

0 

1.3 ± 0.6 

32 Si 

9 

1 

2 

8 ± 3 

31 Si 

12 

3 

4 

18 ± 8 

30 Si 

31 

35 

35 

49 ± 11 

29 Si 

29 

38 

38 

41 ± 10 

28 Si 

21 

26 

25 

21 ± 8 

27 Si 

.4 

.2 

.2 

2.3 ± 1 

32 A1 

0.3 

0 

0 

0.08 ± 0.04 

31 A1 

1.1 

0 

0 

1.1 ± 0.4 

30 A1 

5.7 

1 

1 

4 ± 1.8 

29 A1 

8 

4 

5 

16 ± 6 

28 A1 

19 

19 

19 

31 ± 10 

27 A1 

22 

24 

24 

44 ± 11 

26 A1 

13 

15 

14 

16 ± 6 


TABLE II. TOTAL CHI-SQUARE VALUES FOR ELEMENTS LISTED IN TABLE I 


Element 

produced 

Chi-square value | 

Hyp ergeomet ric 
distribution 

Geometric GDR 
distribution 

Quantum mechanical 
GDR distribution 

Sulfur 

6.1 

14.2 

13.9 

Phosphorus 

12.3 

16.3 

13.3 

Silicon 

13.1 

24.2 

23.7 

Aluminum 

38.6 

21.3 

20.7 
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Figure 1. Production cross sections for isotopes of 
sulfur, phosphorus, silicon, and aluminum pro- 
duced by fragmentations of 213 MeV/nucleon 40 Ar 
projectiles in carbon targets. Theoretical esti- 
mates using the hypergeometric distribution are dis- 
played (dots) and compared with estimates obtained 
from the geometric giant dipole resonance distribu- 
tion (dashed curves) and experimental data from 
Lawrence Berkeley Laboratory (error bars). 
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